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Abstract 

The model sequence calculations by the writers 

of early solar evolution have been extended to the case 

of a varying-G cosmology. In the varying-G sun there 

is less lithium depletion than in a constant-G sun, and 

the time required to burn out central hydrogen is 

decreased by less than the ratio of luminosities owing 

to the much greater total helium production in the 

varying-G sun. At the present time the varying-G model 

8 

has four times the B neutrino flux of the constant-G 
model, a prediction which can be subjected to an experi- 
mental check. - 
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Introduction 
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There have been a number of attempts to incorporate Mach ' s 
Principle into general relativity. Perhaps the most promising 

f 

of these is the attempt by R.H. Dicke and his colleagues to 
introduce a scalar field into general relativity in addition 
to the usual tensor field (Dicke 1964) . Such a theory can be 
formulated in a variety of ways, but the most useful way appears 
to be a formulation in which particle masses remain unchanged 
with time, but in which the gravitational coupling coefficient G 
varies with time. The required rate of variation of G is so slcv 
that the possibility has not yet been subjected to a definitive 
laboratory test. Consequently, it has been necessary to look for 
possible astronomical or geophysical tests of the idea of a varying 


G. 

/ It appears that variations in G will have a great effect cn 

; stellar evolution. G enters into the equations of stellar struct ure 

I i 

^ ' i 

j through the equation of hydrostatic equilibrium: 

I ^ = GM(r)o 

dr ■ ^2 (1) 

j 

i ■■ '■ 

where M(r) is the mass inside radial distance r, p is the density 

I ' and P is the pressure. When 6 is 

' 1 : . I 

) ^ ; ■ ; : 

( ' ^ 

I . . ■ j . ' 
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the pressure gradient is correspondingly larger, and hence any 
structures calculated consistent with the equation of hydrostatic 
equilibrium are more compact. This means that the central tempera- 
tures of the stars would be larger, the opacity would be; consider- 
ably reduced, and stellar luminosities are likely to be much 
greater thaui at the present time. In fact, Dicke has estimated 
that the luminosity of a star in radiative equilibrium might well 
vary as the seventh or eighth power of G. 

Our present evolutionary study was motivated by two possible 

i effects which a varying G would have upon the evolution of the 

I 

I sun. These effects have concerned the rate of destruction of 

■! 

' lithium in the outer layers of the sun and the total' evolutionary 
lifetime of the star of solar mass. 

The essentials of the lithium problem are as follows. Lithium 
is observed to be depleted in the sun, relative to lithium in the 

j meteorites, by at least 2 orders of magnitude (Shima and Honda 1963, 

i Sxiess and Urey 1958) . The question then arises as to %iAiether the 
sun Or the meteorites is normal. Some other stars are observed 
to have much more lithium than the sun; in fact their lithium 
abundances are quite comparable to those in meteorites. Cameron 
(1962) has suggested that the meteoritic abundance of lithium is 
normal, euid that the sun and certain other stars have been depleted 
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in lithium to a great extent. On the other hand. Fowler, Greens tein, 
and Hoyle (1962) have suggested that the normal lithium abundance 
is extremely small in the matter out of which the solar syst.c;ir: 
condensed, and that the amount present in the stn and in the 
meteorites was made by spallation and neutron capture processes 
taking place in the early history of the solar system. Thus it 
is of great interest to determine to what extent lithium would be 
depleted in the normal course of solar evolution by therraonucb ar 
reactions . 

We have recently shown (Ezer and Cameron 1965) that the e- ly 
phases of evolution of the sun involve a fully convective objera 
which contracts toward the main sequence. After some time, the 
central regions of the primordial sun change over to radiative 
equilibrium, and the lower boundary of the convection zone iipp- aches 
the surface. There is still an outer convection zone in the s. , 
indicating that the sun has never been fully radiative. Lithiu::> 
has certainly been destroyed toward the center of the sun, so that 
the question occurs as to whether it would ever be destroyed i; 
that part of the outer layers which constitute the convective zone. 

g 

Our recent study indicated that while the depletion of Li was 
quite large during the contraction phase, the depletion of Li^ 
was relatively minor, amounting to a factor of only 1.5. 



There are many processes by which the depletion of lithium 
in the sun could have been more extensive than indicated in the 
above calculations. These processes include the effects of 
meridional circulation owing to a rapid early rotation of the sun 
and a possible extensive mass loss of the outer layers of the sun 
Meridional circulation would allow interchange of matter between 
the convective envelope and the inner parts of the sun which are 
in radiative equilibrium. Mass loss from the surface of the sun 
would carry away lithium into space; the remaining lithium would 
be diluted as the outer convection zone then extends into deeper 
layers. Neither of these questions have yet been quantitatively 
examined. However, Professor Dicke suggested to us that in the 
early contraction phases of the sun, if G is greater than at 
present, the interior luminosity of the sun would be also greatly 
increased, and this might cause the outer convection zone to exte 
to a greater depth at any given time during the contraction phase 
His suggestion was therefore that a more extensive convection 
region in the early contraction phases would greatly increase 
the lithium depletion. This suggestion was one of the motivating 
factors in our decision to study solar evolution with varying G. 
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The question of the solar evolution lifetime is intimately 

connected with the abundance of helium in the sun and in the 

galaxy as a whole. In an expanding model universe of the typ': 

that Dicke suggests to be consistent with the variation of G, t:: 

9 

present age of the universe is about 7 x 10 years, and hence th 

sun and solar system would have to be formed at the relatively 

9 

early universe age of 2.5 x 10 years. This raises important 
questions concerning nucleosynthesis . If the universe goes 
through a high temperature-high density stage, then there is the 
possibility that the large abundance of helium in the sun and 
stars is of primordial origin (Peebles 1965). The alternate 
possibility that the helium has been manufactured in stellar 
interiors was recently considered by Truran, Hansen, and Gamer or 
(1965). Their conclusion was that inadequate amounts of heli..r.'i 
could be made in massive rapidly evolving stars, but that there 
was a possibility that sufficient amounts of helium might be mac 
in stars of approximately solar mass. Most of the mass of the 
galaxy must be processed through such stars in order to manu- 
facture the initial solar abundance of helium. The crucial 
question regarding the possibility of making helium in a more 
rapidly evolving universe, such as suggested by Dicke, thus 
involves the question of whether the lifetime of the sun would 
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be shortened to much less than 2.5 x 10 if the sun was 
formed rather early in the lifetime of the universe. Professor 
Diclce had suggested that such a shori-eneu lUretiine would occur 
as a result of the increase of G at earlier times. This has been 
a second motivating factor in the investigations reported in this 
paper . 


The Variation of G 


The general relativistic and cosmological equations which 

I 

! incorporate the scalar and tensor gravitational fields have been 
. discussed by Brans and Dicke (1961) . In such a cosmology G varies 
very rapidly in the early stages but much more slowly at later 
times. These variations of G must be determined by numerical 
integrations of equations 48 to 58 of the article of Brans and 
Dicke. During the earlier phases of the expansion an analytical 
solution may be used, which is given by equations 60 to 63 of the 
paper by Brans and Dicke. Professor Dicke has recommended to us 
the following sets of parameters which give a reasonable model of 
the expanding universe consistent with present experimental data: 

cu = 5 

9 

present time = 7.0 x 10 years 

... 2 

decelleration parameter - a a/a = 


1.25 


4 
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which yields at the present time 

G/G = — 2.01 X 10 years ^ 
o 

Tjj = a/a = Hubble age = 13.5 x 10^ years. 
Professor Dicke has provided us with numerical solutions to 
cosmological equations for these parameters which have been 
obtained by Professor Peter Roll. Variation of G with time 
corresponding to this model is shown in Figure 1. 

We have approximated the variation of G by the followij 
expressions: 


r - m 

o 9 


0.1083 


for the interval 


0 < T_ ^ 0.905; 

y 

^ = 1.3943 - 0.1507 T + 0.00664 + 0.00499 

G y y 

o 


for the interval 


0.905 ^ Tg ^ 2.50; 


= 1.25130 - 0.05116 Tg + 0.00218 Tg^ 
o 


for the interval 


Tg s 2.50 


\><o 
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Results and Discussion 

Evolutionary sequences of solar models were constiucted in 
the same manner as described in our previous paper {Ezer and 
Cameron 1965) . The only modification used was the treatment of 
G as a variable rather than as a constant in the equations. 

There is a conceptual error associated with this treatment 
which, however, does not make a significant error in the final 
results . 

The gravitational potential energy of the sun is proportional 
to G. According to the virial theorem, the binding energy of the 
sun is half the gravitational potential energy. Hence, as G 
decreases with time, the binding energy of the sun should also 
decrease with time. The only source of energy to supply this 
decrease in the binding energy is the scalar gravitational field 
itself. Consequently, a complete theory of stellar structure in 
the presence of varying G would include an energy input to the 
sun due to the variation of the scalar field. Such an energy input 
term was not included because the nuclear energy released in the 
solar interior during the sun's lifetime is orders of magnitude 
greater than the scalar field variation energy input. Consequently, 
the neglect of this extra term can have very little effect on 
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the final results. 

The evolution of the sun was started at the universe age 

9 

T- = 2.5 X 10 years. The evoliatxon v.as started at tne 
0 

of stability, at which the gravitational potenti.al energy of . 
sun is just sufficient to supply the internal thermal energy, 
the dissociation energy of hydrogen molecules, and the necessary 
amount of ionization energy of hydrogen and helium in the inter or. 
The solar models were followed through the contraction stage or. o 
the main sequence, and then through the period of hydrogen burn, ng 
until hydrogen was just exhausted at the center. 

The characteristics of representative models, are given in 
Table 1. In this table the first column gives the solar evolution 
time in years measured from the threshold of stability. The 
second and third columns give the radius and luminosity of the 

I 

model in present solar units. The fourth column gives the log ithm 
of the effective temperature of the model, and the fifth columr 
gives the photospheric density in grams per cubic centimeter . the 
remaining columns give, respectively, the central density in gr ns 
per cubic centimeter, the central temperature in °K, the fraccicnal 
mass of hydrogen at the center of the sun, and the relative conrri- 
bution of the gravitational energy release to the total energy 
flowing out of the sun. 
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The luminosity is plotted as a function of .the effective 

temperature in Figure 2. In this figure the dashed line represents 

the track which has previously been found for the sun without 

any variation of G; this track has here been extended to the point 

where the central hydrogen has been depleted in the sun at a time 

of 10^^ years, and through the subsequent formation of a hydrogen- 

10 

burning shell source aund evolution up to an age of 1.4 x 10 
years. The solid line shows the track followed by the sun with 
variation of G. It may be seen that this line lies at slightly 
higher effective temperatures than in the constant G case. Further- 
more, the solid line turns toward the upper left at a higher 
luminosity than for constant G, and it intersects the main sequence 
position also at a higher luminosity. After reaching the main 

sequence, the line curves down toward the lower right, and at an 

9 

age of 4.5 x 10 years it has about the same effectxve temperature 

as the sun, but the luminosity is about 20% too high. Hydrogen 

9 

is exhausted at the center at an age of 9.8 x 10 years; the 
tracks turns up slightly before this point is reached. 

As in the case with constant G, the nearly vertical contrac- 
tion path of the early sun represents models which are fully 
convective. The tracks turn up toward the lower left after a 
central radiative core has formed in the models and as this 
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radia'tiv’e core grows to take in most of the mass of the star . 

After the sun has become mostly radiative, a small central convec- 
tive core forms, grows until a ma>:imum size is reached, and tran 
gradually disappears. The characteristics at the boundaries of 
the radiative and convective regions are shown in Table 2. In 
this table the first column again gives the solar evolution time 
in years. The other columns give the characteristics of the boundary 
between a radiative and a convective region. The first sequence 
of values refers to the top of the central radiative core; the 
later values refer to the top of the central convective core ent* -r 
the bulk of the star has become radiative. The second and subsi...quent 
columns give the fractional mass of the core at the boundary, the 
fractional radius of the boundary, the boundary temperature in " K, 
and the boundary density in grams per cubic centimeter. 

This table may be compared with Table 6a of Ezer and Ccmci . n 
(1965). Such a comparison indicates that, relative to the models 
with constant G, the models with varying G start to develop a 
radiative core at the center sooner. At a given luminosity the 
models have a larger radius and a lower central density. Tne 
resulting low opacities cause a development of a radiative core 
at the center of the sun at considerably lower temperatures than 
for the constant G cases. The temperature at the bottom of the 
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convective zone decreases as the convective zone 'recedes toward 

the surface; the temperature thus never reaches high values. 

This behavior results in less lithium burning in the varying 

G models than in the constant G models. V/e ha^'^e determined by an 

7 

approximate numerical calculation that Li is depleted by only 
10% in the outer convective zone of the sun with varying G. Thus 
our first significant conclusion is that a variety of effects 
associated with the complicated behavior of stellar opacities 
conspire to invalidate the simple prediction which Dicke made 
concerning lithium burning. 

i 

; The decrease of the luminosity and of the effective temperature 
during the main sequence hydrogen burning phase of the sun with 
varying G is an effect which has already been found by Pochoda and 
Schwarzschild (1964). It may be noted that the models have too 

9 

large a luminosity by some 20% at the sun's age of 4.5 x 10 years. 
We estimate that the luminosity would have come out correctly if 
the helium content by mass of the initial sun had been reduced 
from the assumed value of 24% to about 22%. This alternate choice 
of helium content would be quite consistent with the uncertainties 
in the hydrogen-to-helium ratio in the sun. 

It is at first sight surprising to find that the models 
with constant G and with varying G both take very nearly 10^^ years 
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to biirn out their central hydrogen, despite the fact that tha models 
with varying G have a much higher luminosity throughout mosc of 
the solar lifetime. This can occur only if much more hydrogen 
has been consumed in the interior of the sun at the time of :e ural 
hydrogen burnout for the varying G case than for the cons La:; u 
case. 

We have investigated this point in several ways. Figure : 

shows the ratios of the net rates of destruction of hydrogen by 

the two branches of the proton-proton chain, r., is the r<;;r . of 
3 7 

the He (a,Y)hi reaction. The branch of the proton-proton cha n 

... 4 

containing this reaction makes one He nucleus per one basic 

proton-proton reaction. The rate r is the rate of the 
3 3 4 

He (He ,2p)He reaction. In this branch of the proton-proton 

4 

chain the formation of one He nucleus requires two basic pro .- 
proton reactions. The ratio r ^/(r_ . + 2r ) increases .viv; 

temperature and He^ content. It may be seen that when this rai lo 
is high, a greater fraction of the luminosity of the sun will he 
supplied by the material at the very center, and hence the ma:. rial 
at the center will burn out relatively sooner. Figure 3 show: 
this ratio for the varying G and constant G cases very early in 
solar evolution on the main sequence. It may be seen that in Lae 
constant G case this ratio is rather small but has a fairly sharp 
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peak toward the center. In the varying G case, the ratio is higher 

and has a much broader flat peak at the center associated with 

larger central convection zone in the sun. 

Figure 4 shows the fraction by mass of helium in the; present 

9 

sun, corresponding to the age 4.5 x 10 years, for both the varying 
G and the constant G cases. It may be seen that much more helium 
has been formed for the varying G case than for the constant G 
case, as was expected. It may also be seen that the excess helium 
extends very significantly to larger mass fractions. Figure 5 
: shows a similar plot for the two cases with the central hydrogen 
just depleted. Here it may be seen that in the constant G case 
the central hydrogen is depleted just at the very center of the star 
and the hydrogen content then rapidly increases away from the 
center. In the varying G cases, however, hydrogen has burned out 
nearly uniformly over a fairly wide region, and hydrogen then 

increases away from this flat plateau. The amount of helium formed 

■: ; 

in the sun with varying G is much greater than that with constant 
G, and this indicates why the two model sequences can take nearly 
the same time to burn out the central hydrogen. 

We have found some significant differences in the spectrum 
of neutrinos emitted from the present sun in the constant G and 
veurying G cases. Table 3 shows the neutrino fluxes at one 
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astronomical unit for the solar models with constant G which were 
reported in our previous paper (Ezer and Cameron 1965). The 
first four columns give the evolutionary age of tlie sun in years, 
the radius in solar units, the luminosity in solar units, and .Ir 

central fractional mass of hydrogen. The subsequent columns give 

1 +2 7 — 7 

the neutrino flux from the reactions H (p,3 v)D , Be (e ,v)Li , 

B®(p%)Be®*, and The last column 

gives the total neutrino flux. Table 4 gives a similar set of 

numbers for the solar models with varying G. In order to compare 

9 

Tables 3 and 4 at 4.5 x 10 years, it would appear appropriate 

to decrease the neutrino fluxes in Table 4 by about 20%. 

On this basis the present sun with varying G emits approx- 

8 

imately four times as many neutrinos from the decay of B as 
does the present sun with constant G. This difference has pare: ilar 

significance in view of the fact that the large solar neutrino 
detection experiment of R. Davis (1964), consisting of a large 
underground tank containing a chlorine compound in which solar 
neutrinos produce radioactive A , is most sensitive to the b' 
neutrinos from the sun. His experiment, which is now in prepara' ion, 
is expected to detect the solar neutrino flux, and therefore thei 
is a real possibility that it will constitute a test of the varia ion 
or nonvariation of G if solar models can be sufficiently refined 
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so that there is agreement about the expected neutrino flux in 

either of the cases. It may be noted that the best prediction 
8 

of the B neutrinos for the present sun with constant G due to 
Sears (1964) gives a value which is interirediate between the two 
numbers given in Tables 3 and 4. 

The failure of the solar model with varying G to burn out 
its central hydrogen significantly faster than the solar model 
with constant G raised questions in our mind concerning the 
possible role of stars of about solar mass in the earlier history 
of the universe in the cosmology in which G varies. Consequently, 

j 

calculations of the evolutionary tracks where repeated, starting 

the sun at the threshold of stability at universe ages 
9 9 

Tq = 0.5 X 10 years and T^ * 0.1 x 10 years. The results are 
displayed in Tables 5 through 8. Tables 5 and 6 are equivalent 
to Tables 1 and 2, respectively, for the models starting with 

9 

Tq = 0.5 X 10 years. Similarly, Tables 7 and 8 are equivalent 

9 

to Tables 1 and 2 for the models starting with T^ = 0.1 x 10 years. 
In these models G is initially larger than in the previously 
calculated models, and consequently the central temperatures in 

9 

the initial sun are higher. For the case T^ = 0.1 x 10 years, 
the carbon-nitrogen-oxygen cycles of nuclear reactions make 
significant contributions to the energy generation in the center. 



The evolutionary tracks and the luminosity-effective te npora- 


ture diagram are shown in Figure 6. It may be seen that the 

models starting at smaller ^universe stjes go to siignificantly 

higher luminosities and effective temperatures than would be tie 

case for the present sun with varying G. The tracks then exhibit 

much more extensive movement toward lower luminosities and effective 

temperatures. The maximum mass fraction of the central conveceve 

9 

core extends to 17% and 28% of the mass for = 0.5 x 10 vc-ar\ 

0 

9 

and 0.1 X 10 years, as compared to 15% of the mass for tne ca - 

9 

Tq =s 2.5 X 10 years. The times required for depletion of t;:e 

9 

central hydrogen in the cases T^^ = 0.5 x 10 years and 0.1 x lu 

9 9 

years are 3.5 x 10 years and 2.4 x 10 years. It may be seen .hat 

the first of these cases has a marked reduction in this burn-or 

9 

time as compared to T^ = 2.5 x 10 years. However, there i n 

> 

much additional reduction in lifetime for the case T^ = 0.1 x 
years, presumably owing to the much greater extent of the cc.ee: 1 

convective core in these models. The more extensive central 
hydrogen burning which takes place in these other model sequenc 
is indicated in Figure 7. It may be seen that the earlier the 
sun starts to evolve in the age of the universe, the greater is 
the burn-out of hydrogen in the interior at the time hydrogen 
becomes exhausted at the center. 



These results are particularly significant in connection with 
a question of helium formation in the galaxy. ?or the ca£-e of 
stellar evolution with constant G, Truran, Hansen, and earner on 
(1965) have shown that massive stars can contribute very little 
helium to the interstellar medium, but there is a chance that the 
large amount of helium in the sun could have been made as a result 
of stellar evolution in stars of approximately solar mass. However, 
the associated universe ages must be rather long and the fraction 
of initial galactic gas that must be processed is large. The 
present results indicate that such stellar helium production is 
very difficult in this varying G cosmology. Thus the large amount 
of helium in the initial sun may have to be attributed to primordial 
processes (Peebles 1965). 
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Evolutionary sequence of solar models avith varying G and beginning 
of evolution at a universe age T = 2.5 x 10® years. 
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TABLE 1 (continued) 
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TABLE 1 (continued 
























TABLE 2 (continued 














Neutrino fluxes at 1 astronomical unit cm sec 
for evolutionary solar models with constant G 
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Evolutionary sequence of solar models with varying G and beginning 
of evolution at a universe age T « 0.1 x 10® years. 
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Figure Captions 
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Figure 1: The variation of the gravitational coupling coefficient G 

”8 2 2 

with time. G(now) = 6.67 x 10 dynes cm /gm . 

Figure 2; i Evolutionary track for the sun. The solid line is 
calculated with G varying in time. The dashed line 

s ■ ■ . , 1 . ' 

i j represents the evolution of the sun with constant G. 

; J . * -f ; 

The crosses on the tracks correspond to the points at 
which the central hydrogen has been depleted in the sun. 
Figtire 3: The ratios of the net rates of destruction of hydrogen 

by two branches of the proton-proton chain, for varying G 
and constant G, very early in the main-sequence solar 

' ^ d ' i • • • ■ ' 

evolution. 

Figure 4: Helium concentration by mass as a function of mass fraction, 
in the present sun, for both varying G and constant G 
cases. 

Helium concentration by mass as a function of mass 
■'- fraction, for varying G and constant G models having 
; central hydrogen just depleted. 

Figure 6: The Ivuninosity-eff active temperature diagram for a 

one solar mass steur starting evolution at the indicated 
■ ? ' universe ages . 




1 Figure 5: 
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Figure 7: The fraction by mass of helium in the sun when 

central hydrogen is just depleted. Each curve 
represents the indicated universe ages for the 
beginning of stellar evolution. 
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Figure 2 
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